MARKOVIAN SOLUTIONS OF INVISCID BURGERS EQUATION 
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Abstract. For solutions of (inviscid, forceless, one dimensional) Burgers equa- 
tion with random initial condition, it is heuristically shown that a stationary 
Feller-Markov property (with respect to the space variable) at some time is 
conserved at later times, and an evolution equation is derived for the infini- 
tesimal generator. Previously known explicit solutions such as Frachebourg- 
Martin's (white noise initial velocity) and Carraro-Duchon's Levy process 
intrinsic-statistical solutions (including Brownian initial velocity) are recov- 
ered as special cases. 



1. Introduction 

We consider the inviscid Burgers equation d t u + d x (^u 2 ) = with random 
initial data uq. Burgers equation has originally been introduced [2] as a ID model 
of turbulence. Although it is now clear that it does not exhibit lots of features of 
"true" turbulence, we nevertheless still think it is a good equation on which one can 
try and find new methods to apply on Euler equation. Having this in mind, taking 
random initial data seems quite a natural problem. It is also physically relevant in 
the contexts of interface dynamics, of aggregation of particles [Jj], and some others. 
Burgers equation with a random force on the r.h.s. has also been studied, mainly 
as a "benchmark" to test methods designed for (Navier-Stokes) forced turbulence, 
many of which turn out to produce spurious predictions when applied to the simpler 
Burgers case. See [5] for instance. 

The case of a Brownian initial data has already been investigated by Sinai 
Carraro and Duchon 0] showed that Levy processes are conserved by Burgers 
equation. They also obtained the explicit evolution equation for the characteristic 
function of the Levy process solutions of Burgers. A noticeable point is that they 
made no use of the Hopf-Cole construction of the solution (Bertoin pQ recovered 
essentially the same result with Hopf-Cole). 

We will essentially follow Carraro and Duchon : first we define what we call a 
statistical solution of Burgers equation, and write an infinite set of equations for the 
71-point functions of such solutions. We show that the assumption that the process 
is Feller (in space) for all time yields an evolution equation for the infinitesimal 
generator of this process. Conversely, a Feller process whose generator satisfies this 
equation is a statistical solution of Burgers equation. This will allow us to recover 
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Carraro and Duchon's result on Levy processes, as a special case. Frachebourg and 
Martin's explicit solution in the case of an initial white noise velocity is also a 
particular solution to our equation. 



2. Notations and definitions 

A Markov process u(x) X £r can be characterized by its one point and its transi- 
tion probabilities p x (du) and q XiV (u,dv) , x < y, that satisfy, Vxo < ■ ■ ■ < x k and 
/, borelian positive, E[n£= fi{ u { x i))] = 

J Px (du )fo{u ) / 9x0,^1 (wo,<fui)/i(ui) ■ • • J q Xk . 1 ,x k {u k - 1 ,du k )f k {u k ). 

A Markov process is homogeneous if its transition probabilities q XzV depend on x 
and y only through y — x. In this case, we write qu instead of q XlX +h ■ 

A process u(x) X £R is stationary if and only if it is translation invariant : the law 
of (u x+Xl , . . . , u x + Xrl ) does not depend on x. Hence a Markov process is stationary if 
and only if it is homogeneous and its one point probability p x (du) does not depend 
on x. 

If u is a homogeneous Markov process, h > 0, and / is a continuous function 
vanishing at infinity, we put Qhf(u) = J f(v)qh(u,dv) . 

A Feller process is a homogeneous Markov process such that for each /, for 
each h > 0, Qhf is also continuous and vanishes at infinity, and lim^o Qhf = f 
pointwise. 

A Feller process always has a cadlag version 

One can define the infinitesimal generator of a Feller process : it is the operator 
A, defined for all the functions / such that the limit below exists, by 

VxeR,Af(x)= lim Qhfix) - f{x) 

Formally, Qh — expQiA) , Q' h := dQ/dh = AQh and an invariant measure po 
satisfies Ap = . 
The 



3. Statistical solutions of Burgers equation 

We will closely follow |3] (see also JHl)- Let E be the space of cadlag real 
functions equipped with the smallest er-algebra C(E) such that for each x E R, 
u i— > u(x) is measurable. Let V be the set of real C°° functions with compact 
support. A probability fi on E is then characterized by its characteristic function 

uEPh / exp [i / u(x)v(x)dx] dfi(u) = jl(v). 
Je Jr 

Let Mo : (il,A, P) — > E be a random process, defined on some probability space, 
and let fio : C(E) — > [0,1] denote its probability law. Assume u(x,t) is a (weak) 
solution of Burgers equation with u(., 0) = uq, u(., t) G E for t > 0, and everything 
makes sense in the following calculation : integrability, and differentiability with 
respect to t. Let fi t denote the law of u(.,t). Formally, one then gets for each 
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v e V : 

dtfit(v) = I d t {exp[i I u{x)v(x)dx]} djj, t {u) 



<9t{exp [i / u(x,t)v(x)dx]} d/j,o(uo) 
t; lit 

exp / u(x,t)v(x)dx)]dt[i / u(x, t)v(x)dx] dfj,o(uo) 
■In In 

exp [i I u(x,t)v(x)dx] i I —u(x,t) 2 v'(x)dxdfio(uo) 
■In In - 

1 ' S2.J, 



= i / —u(x) v {x)dx exp [i / uv]dfj,t(u) 

J E J R 2 </ R 

Hence our definition of a statistical solution of Burgers equation : 

Definition 1. A statistical solution of Burgers equation is a set (nt)t>o of proba- 
bilities on {E,C{E)) such that for any v £ T>, 

(i) dtfit{v) = i I I -u(x) 2 v'(x)dxexp[i j uv]dfi t {u) 

Je Jr 2 J R 



Let us assume now that we have a statistical solution of Burgers equation, 
(A t t)t>o , and that for all t, all the moments of fit are well defined. Then one 
can write 

exp [i / u(x)v(x)dx] = — / J^J u(xj)v(xj)dxj 
Equation |QJ thus becomes Vw 6 D, 

oo . n ~ n 

(2) l / a t E[n u (^x^)] n d ^ = 

= 2 — / E[]^ w(xj)z;(xj)ii(x) 2 t' / (x)]dx = 

00 „ n w \ 71 

= ^ ! " E E M^)^Hv II «(*fc)«0»*)] 11 dx k 

4. Evolution equation for Markov solutions 

We are now looking for solutions such that at each time t, x 1— > u(x,t) is a 
stationary Feller process (with respect to space x). We are going to show that 
for such processes, the infinite set of equations @ is equivalent to an evolution 
equation for the infinitesimal generator of u. 

We thus assume now that the solution x 1— > u(x, t) is a stationary Feller process, 
with one point probability p(du,t) and transition probability qh{u\,du2,t), the 
equation J2J) becomes Vu £ T> : 
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oo „ 

(3) 



dx\ . . . dx n u\v(xi) . . . u n v(x n ) 

< — <x„ 



x dt[p{du 1 )q h2 (ui,du 2 ) ■■■qh n (u n -i, du n )\ 

oo „ 

y^i n dx ... dx n p(du )q hl (it , dm)... q hn (u n -i,du n ) 

i J x n <-<x„ 



V V 

x [u — (x ) H h u„ — (x„) u v(x ) ■ ■ ■ u n v(x n ) 

v v 



oo „ 

= * X/™ / dx ... dx n u v(x a ) . . . u n v(x n ) 

n=l Jx <-<x n 

x p(du )q hl (u ,dui) ...qh n {u n -i,du n ) 
x {ui^2-(ui, du 2 ) + u 2 [^(u 2l du 3 ) - ^{u u du 2 )} + ... 

Qh 2 qh 3 qh 2 

+ u„_i[— (u„_i,(iu„) - />n ~ 1 {u n - 2 ,du n -i)] - u„^.(u„_i,rfu„)} 

qh n qhn-i qh n 

(by integrating by parts; we note hj — Xj — Xj-i and q' h — dqh/dh). 

This equality is equivalent to the following infinite set of equations : Vn e 
N*,Vxi < ... <x n : 



(4) 2d t E[u(x 1 )...u{x n )} = 

u x . . .u n p{du 1 )q h2 (u 1 ,du 2 ) . . . q hn {u n -i, du n ) 



/ 



x {ui — (mi, du 2 ) + u 2 [—(u 2 ,du 3 ) - —(u l7 du 2 )] + . . . 

qh 2 qh 3 qh 2 

+ u n -i[^ hL (u n -i,du n ) - <lh "~ 1 (u n - 2 ,du n -i)] - u n — (w„_i,rfu„)} 

qh„ qhn-i qhn 

One then gets the evolution equations for p, q and A by taking limits in which 
some of the x^s are equal. If one makes every Xi tend to Xi, the preceding set of 
equations gives formally, VneN* : 

2 J d tP {du)u n = J p(du)(-UAU n + U n AU){u) 

where U n denotes the function iihii". 

If one makes some of the x^s tend to xi, and the others tend to x 2 = x\ + h, one 
then gets Vn e N*,Vfc < n,Vxi £R,Vli£ R+* : 

(5) 2d t E[u(x 1 ) k u(x 1 + h) n - k ] = J p(du){-UA(U k Q h U n - k ) + U k+1 Q h AU n - k 
+ U k [A(UQ h U n - k ) - AQ h U n - k+1 + Q h {U n - k AU) - Q h (UAU n - k )]}(u) 
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One then easily finds, if r\ is in the domain of A : 

(6) 2 J d t p(du)r](u) = J ' p(du)[-u Ar](u) + rj(u) AU (u)] 

(7) 2d t QhV = UAQhV + A(UQ h f})-Qh(UAri) 

-AQ h {Urj) + Q h (vAU) - AUQ hV 
These two equalities sum up into one : Vry in the domain of A, 

(8) 2d t Ar) = U A 2 tj - A 2 (U rj) + A(r) AU) - AU Arj 

or, introducing the operators Mjj and M A u defined as Mur){u) = ur)(u) and 
M AU r](u) = AU{u) V {u) : 

(9) 2d t A = M V A 2 - A 2 M V + AM AU - M AU A 

If this latter equality holds, one can easily check that if Ap = for all time, then 
p verifies ©, and Qh = exp(hA) verifies @. 

Hence a Feller statistical solution of is solution of |JBJ and @, which are 
equivalent to 0. 

Conversely, it is a matter of simple algebra to check that © and 10 imply J2J : 
indeed one can then write for any x\ < • • • < x n (recall hi = x% — Xi—x) : 



(10) 2d t E[u{xx) . . . u(x n )} = 2d t J P (du)MuQ h2 . . . M v Q hn U{u) 
= 2 [ dtp{du)MuQ h2 . . . MuQ hn U(u) 



n - 

2 H I P( du "> M uQh 2 ■ ■ ■ MuQhj-iMu d t Q h] MuQ h]+l . . . MuQ hn U(u) 
P (du)u[-AMuQh 2 Mu . . . Q hn U(u) + Q h2 M v . . . Q hn U(u) AU{u)] 



n „ 

+ / P( du ^ M uQh 2 ■ ■ ■ MulMuAQh^j + A(U Q hj nj) - Q hj (U Aty) 

- AQ hj (U + Q hj ( Vj AU) - AU A m ] 
where rjj = MuQh j+1 ■ ■ ■ MuQh n U. Many terms cancel, one gets 

= / P (du)uJ2Q h2 Mu . . . Qh^MuIMuQ'k. - (/,, M, '/., 

which is just one integration by parts away from J3J|. 

Therefore, if u(x,t) is a Feller process, it is a statistical solution of Burgers if 
and only if its infinitesimal generator is solution of JjJ}. In some sense, the Feller 
assumption yields an exact closure of the infinite set Of course, nothing guar- 
antees the existence of solutions of @, although we show later that the Brownian 
and white noise initial cases give formal solutions to it. Nevertheless, a close look 
at Bertoin's proof using Hopf-Cole pQ makes us strongly suspect that the absence 
of positive jumps may be essential to guarantee the existence of solutions. This 
would also be reasonable from a physical point of view : solutions with positive 
jumps are unphysical. 
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5. The case of Levy processes 

We will see how one can recover formally the results of 0|. The initial velocity 
ito is here supposed to be a Levy process (which means that it has independent and 
stationary increments) of finite variance having no negative jumps. This covers 
in particular the case of uq Brownian. Such processes are characterized by their 
second exponent <j>, defined by Vie < y,VA <E R + : 

£{exp [X(u (y) - w (x))]} = exp[(y - x)(f>(X)} 

A Levy process can be considered as a limit case of stationary Markov process (the 
one point distribution p is replaced with Lebesgue measure). One can also formally 
define an infinitesimal generator by the relations : V A € R + , 

where we have noted e\ the function u t— > exp(Au) (which of course is not in the 
domain of A . . .). One can inject these relations into the evolution equation JJjJ. 
Using ue\{u) = d\e\(u), and AU = constant, one gets an evolution equation for 
4> ; it turns out that this equation is also the Burgers equation : 

(11) 2d t <f>(X) = -d x 4> 2 

Carraro and Duchon ^ have checked that if 4>o is the exponent of a Levy process 
of finite variance with negative jumps, (jf l|l has a smooth solution for all time t > 0, 
which is still the exponent of a homogeneous Levy process with negative jumps. 

Hence such Levy processes are conserved by the Burgers equation. The Brownian 
case corresponds to (/>o(X) = A 2 /2, and this yields <j){\ t) = (1 + At — \/l + 2Xt)/t 2 . 

6. Evolution equation for the jump process 

The infinitesimal generator of an arbitrary Markov process can be written as the 
sum of three terms (see [H]) : a diffusion term, a drift term and a jump term : 

Af(x) = a(x)f"(x) + b(x)f'(x) + J n(x, dy)(f(y) - f{x)) 

The measure n(x, dy) represents the jump part of the process : it gives the number 
of jumps going from x to y. In our case, all these coefficients will of course depend 
on time. To write an evolution equation for n , we assume b = 1/t and a = for 
t > , and all jumps are negative. Equation Q then yields an evolution equation 
for n : Vx > y, 

(12) 2d t n(x, dy, t) — j(x - y)(d x n(x, dy, t) - d y n{x, dy, t)) 

n(x, du, t) [(x — y)n{u, dy, t) + (y — u)n(x, dy, t)] 

y 

(x — u)n(y, du, t)n(x, dy, t) 

-oo 

7. The case of an initial white noise process 

Frachebourg and Martin |7] have investigated the case of an initial white noise ve- 
locity. Using the Hopf-Cole construction, they obtain explicit formulas for the laws 
of u(x, t) and its jumps. They actually rederived results about Brownian motion 
with a parabolic drift that had been previously established by Groeneboom |H] out 
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of the Burgers context. Using Frachebourg and Martin's results or Groeneboom's 
paper, the infinitesimal generator in the case of an initial white noise process is 
found to be, in the case where < uo(x)uo(y) >= (X/8)S(x — y) : 

Af(x) = -f{x) + 4 J Jf(y) - f(x))(x - y)-^!^ 3 - yt^)dy 

where / and J are given by their Fourier and Laplace transforms in terms of the 
Airy function Ai : 

-r, \ 1 P°° , exp(uz) 

(13) J(u) = — / dz- 



2iTT J^ loc 2 1 /3Ai(2- 1 /3 z ) 

1 r ioc 9 2 / 3 AiV9-V3^ 

,i4) 2/w - ^-*+±L^* +l2 > m 

We have checked that the evolution equation 10 is indeed verified : it amounts 
to expressing convolutions like ul * J , ul * ul, ul * u J in terms of J 1 and It can 
be done using relations (fH5|l and (fTljl and the fact that Ai' (a:) = xAi[x). 



8. Conclusion 

We have heuristically shown that for Feller stationary processes, Burgers equa- 
tion is equivalent to an evolution equation for their infinitesimal generators. It gives 
strong evidence that the Feller property is conserved by Burgers equation, although 
we suspect that the negativity of jumps in the initial velocity should be required. 
Our evolution equation provides an equation for the jump process, and this might 
lead to other exact statistical solutions of Burgers equation. The previous exact 
solutions concerning an initial Brownian or white noise velocity are both particular 
solutions of our equation. 
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